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Introduction t< 
Redirect Missi 


■ The Asteroid Robotic Redirect Mission (ARRM) is a mission concept to 
retrieve a large sample of asteroidal material, either: 

- Option A: a small Near-Earth Asteroid (NEA) via inflatable bag 

- Option B: or a piece of a large NEA via articulated grippers 

■ ARRM will also test large-scale solar arrays and electric propulsion 
systems 

■ ARRM is intended to be a stepping stone to the very large propulsion 
and power systems that will be needed for human exploration of the 
solar system 
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Why ARRM 


■ Mars sample return has been a tantalizing but unattainable goal since 
the beginning of planetary exploration 

■ A very interesting precursor mission might be to collect samples from 
Deimos or Phobos 

- These bodies are windows into both Mars itself and into the primitive 
state of the solar system 

■ Previous sample return studies have focused on retrieving small 
amounts of material, but what if a large multi-ton sample could be 
returned? 

■ The ARRM Option B spacecraft, with its articulated grippers and landing 
legs, could also be used to retrieve a large sample from Deimos or 
Phobos 

■ Alternatively, the same spacecraft could be used to retrieve a sample 
container launched from the surface of Mars 
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Guiding Qu 


■ “How much mass can we retrieve from Deimos, Phobos, or low-Mars 
orbit (LMO)?” 
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Mission Ar 


■ Launch from Earth on either Falcon Heavy (FH) or Space Launch 
System (SLS) 

- FH launch is followed by a double lunar flyby maneuver which results 
in ejection from the Earth-Moon system with a C3 of 2 km 2 /s 2 

- SLS directly delivers the spacecraft to a heliocentric trajectory with a 
C3 of up to 27 km 2 /s 2 

■ Spacecraft travels to Mars and arrives at the sphere of influence (SOI) 
with C3 of 0 km 2 /s 2 

■ Spacecraft descends to the orbit of interest (Deimos, Phobos, LMO) 

■ Spacecraft acquires sample (stay time of 200-600 days) 

■ Spacecraft ascends to the Martian SOI 

■ Spacecraft returns to Earth-Moon system with arrival C3 of 2 km 2 /s 2 

■ Earth flybys on the way to or from Mars are considered 
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Trajectory D 


■ 


Power Assumptions 

BOL Power at 1 AU 

51 kW 

Array performance model 

TJGA 

Array decay rate 

1% per year 

Spacecraft bus power 

0.5 kW 

Power Margin 

15% 

Propulsion Assumptions 

Thruster model 

Fixed I ap and system efficiency 

Input power bounds 

0-40 kW 

Thruster performance 

2000 s 7, S p with 55% efficiency 

Duty cycle 

3000 s I ap with 60% efficiency 
90% 

Propellant margin 

11 % 

Launch Assumptions 

Maximum departure hyperbolic excess energy (C3) 

2 km 2 /s 2 (FH) 

Maximum launch mass 

27 fern 2 /* 2 (SLS) 

13153 kg (FH with 2000 s I, p ) 

Declination of Launch Asymptote (DLA) bounds 

13169 kg (FH with 3000 s I sp ) 
16072 kg (SLS) 

+/ - 28.5° 

Post-launch forced coast 

30 days (SLS) 

Mars stay assumptions 

Desired altitude 

none (FH) 

5844 Arm (Phobos) 

Stay time at desired altitude 

20003 km (Deimos) 

400 (LMO) 

200-600 days 

Earth Return Assumptions 

Maximum arrival C3 

2.0 Arm 2 /s 2 

Post-mission An 

100 m/s 

Minimum spacecraft final mass 

5170 kg (2000 s I sp case) 


5020 kg (3000 s I sp case) 


6 # 
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Full-Mission 
Return Probl 


■ When we first began this project: 

- We modeled the Earth-Mars and Mars-Earth phases of the mission separately 

- We varied the pickup mass and stay time at Mars parametrically 

- We computed the descent and ascent spiral flight time and propellant mass in 
post-processing 

■ But is there a better way? The real question is how much mass can we bring 
home - can the optimizer operate on pickup mass directly? 

■ We constructed a full sample-return model which captures the relationship 
between pickup mass and flight time 

■ The combination of variable pickup mass and fully-modeled descent and ascent 
spirals is highly complex - “time is mass and mass is time” 

■ The full-mission sample return optimization framework presented here is 
implemented in the Evolutionary Mission Trajectory Generator (EMTG), 
Goddard’s open-source interplanetary trajectory optimization tool 
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Multiple Gra 1 
(MGALT) via 


■ Break mission into phases. Each phase starts and ends at a body. 

■ Sims-Flanagan Transcription 

- Break phases into time steps 

- Insert a small impulse in the center of each 
time step, with bounded magnitude 

- Optimizer Chooses: 

■ Launch date 

- For each phase: 

■ Initial velocity vector 

■ Flight time 

■ Thrust-impulse vector at each time step 

■ Mass at the end of the phase 

■ Terminal velocity vector 

■ Assume two-body force model; propagate by solving Kepler’s problem 

■ Propagate forward and backward from phase endpoints to a “match point” 

■ Enforce nonlinear state continuity constraints at match point 

■ Enforce nonlinear velocity magnitude and altitude constraints at flyby 
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Spiral Trajectc 
Edelbaum s M 


■ Assume: 

- The initial and final orbits are circular and coplanar 

- The thrust available to the spacecraft is relatively low 

- The thruster is operated continuously during the spiral orbit transfer and thrust 
is always tangent to the velocity vector 


A I planet 
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■ Note: For ascent spirals, m be f 0re _ spira i is dependent on any mass added at the 
bottom of the gravity well, hence “time is mass and mass is time” 
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0 50 100 150 200 250 


Decision Variable 

This is the standard Sims-Flanagan Jacobian sparsity pattern for a two-phase 
mission with bounded flight time, a bounded return date, and a dry mass constraint 
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Sample Ret 
Low-Thrust 



and 



■ Dependency of flight time constraint on the launch date (launch date 
influences arrival at Mars, which influences spiral time, which influences 
flight time) 

■ Dependency of Earth return date constraint on outbound flight duration 


NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 

ft a.i. solutions 















Sample Ret 
and Low-Thr 



Dependency of return-phase match point constraints on pickup mass 

Dependency of Earth return date and total flight time on pickup mass 

Dependency of dry mass constraint on pickup mass due to its affect on 
propellant used for the ascent spiral 
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Nonlinear Pr 




Minimize f(x) 
Subject to: 

%lb — % — %ub 

c(x ) < 0 

Ax < 0 


where: 

x w , x ub are lower and upper bounds on the decision variables 

c(x ) is a vector of nonlinear constraints 
Ax is a vector of linear constraints 


We use SNOPT as our NLP problem solver 
But all NLP solvers require an initial guess.. 
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Monotonic 
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Results 
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An interestin 
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Another interes 
Deimos Sample 


th flyby 


le8 


x (km) 


-3-2-10 1 2 


departure 


begin spiral 

Mars 


Mars 

5/22/2023 


9/30/2022 

C3 = 0.00 km * * /s J 


Av = 1.08 km/s 

Av = 0.00 km/s 


m = 16171 kg 

..... 

m = 15589 kg 

t-;—£— m 




LT rendezvous 
Mars 

10/8/2021 
m = 8626 kg 


unpowered flyby 
Earth 

6/25/2020 
v x — 2.77 km/s 
DEC = 20.5 
altitude = 7038 km 
m = 9792 kg 


launch 

Earth 

9/28/2019 
C3 = 2.00 km 2 /s' 
DLA = 0.0° 
m = 10244 kg 


intercept 
Earth 
9/20/2025 
v x — 1.41 km/s 
DEC = -30.0 
m = 13439 kg 


end spiral 
Mars 

3/14/2022 
Av = 1.08 km/s 
m = 8316 kg 


-1 


-2 



— 2 , 


NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 

(•) a.i.solutions 


17 * 


le8 y (km) 























Conclusions 


■ We constructed a full-mission optimization framework for sample return 
problems which includes both variable pickup mass and low thrust 
descent/ascent spiral modeling 

■ Our approach allows pickup mass to be optimized directly and captures the 
dependencies of all of the problem constraints on the pickup mass 

■ This technique is efficient and reduces the human work-load necessary to 
design sample return missions without requiring large parametric runs 

■ The complexity of the problem increases, but this complexity is well-handled by 
the autonomous optimization algorithms in EMTG 

■ Using this technique, we found viable applications of the ARRM spacecraft to 
Deimos, Phobos, and Mars surface sample return in the 2020s 
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EMTG is available open-source at 
https://sourceforae.net/proiects/enntq/ 
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